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A method is suggested for predicting the stability of auto- 
mat ‘ caJ.Iy controlled aircraft by a comparison, of calculated f re queue; 
resncnce curves for the aircraft and exporirirn bally determined 
j. .’eqvency -response curves for the automatic pilot . The method is 
applied only to stabilization in roll. The method is expected to be 
useful as a means of establishing the specifications of the perform- 
ance* required of the automatic control, device for pilotless aircraft 
designed as missiles. 


IHTRODCJCTION 


Experience Lao shown that the provision of automatic stabi . 


zation for snail pilotlesc 


urc rax ■ 


aei 


rlgned as missiles is extremely 


difficult. The di:l 


[cult 


is a result of the high-frequency oecj 


la t ions of snail- size aircraft- that require rapid control movements 
an i small time lags, characteristics which are difficult to obtain, 
■particularly when the space available for the control servomotors 
end intelligence units is considered. In an unpublished analysis 


ale at the Langley Memorial Aeronautical Laboratory ox 


MCA, 


the problem of determining the stability of on automatically controlled 
aircraft with lag in the control system was analyzed theoretically by 
k : osuming a simplified equation for the control motion, this equation, 
being obtained from the knowledge of the behavior of the automatic 
pilot. Because of the irregular response char act or i sties often found 
in automtic pilots, however, the control motion is difficult to 
represent mathemati ca.lly and , hence , the simplifies equations of cue 
control wore found to be inadequate for the analysis. 


I'he present paper surge i 


;ts a method for predicting 


stability 


of p.v aircraft bused on the experimental determination of the 
characteristics of its automatic pilot. The procedure consists 


EEG , j T RIC r .CED 



control motion re 
of : in it amplitude 
freedom being inspected. The motion of the cont 
a range of frequencies; the phase c-ngie cf ore c 
the ratio of ampiitin.o of control motion to sir; 
plotted as a function of fiornicney. Similar rnr 
for the autopilot 07 oscillating it and '.-‘ecoixUn 
The two sets of do/or .-.re then compared to dotovr 
airplane will be sf'-blc under control of tee nu J 
method is developed in detail only fox* stabilise 
nay be used by the uirpi'ine designer for cither 
suitability of cr existing automatic pilot for a 
cation or specifying die eh;n*acteristieo cf the 
needed for the application. 


essentially in calc- -lab ring the 
a continuous sinusoidal motion 


O 


T1S0I3 


ra macs of aii-plane, sings 

kj- radius of gyration of airplane r lc v i; Ion pit. 

q dynar-dc pressure, pounds per ecu;.' re foot 

5 wing area , 0 quo re feet 

b wing span, feet 

C ^ rolling- moment coefficient (Tolling none rut 

0 angle of bank, radians 

P angular velocity in bank, radians per cecor; 

6 deflection of aileron, radians 

C ^ rate of change of rolling- moment coefficien 
P velocity in bank, per radian (ft 7 /or) 

C ? rate of change of C 7 with u, nor radian 
o 6 

D differential operator (d/dt) 

o angular frequency, Indians pe?- second 

0 phuso angle (positive value means lead of 
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quj. rod to maintain 
for the degree of 
rod is obtained for 
orh rol motion and 
lane motion are 
■70s are established 
7 die control motion. 
I: if; whether the 
omr ti c pi lot . The 
.ti.cn in roll. It 
lei. or: lining the 
Particular appli- 
a a t erratic pilot 


udinal axis, feet 


71) 


a (dp/at) 


.3 with rngiilar 

cm /is) 


qIiop .4 of 0 
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ma:d.num amplitude of 0 


l 


M- 



t(t) 


control-amplitude ratio (ratio of control deflection to 
airplane di splucenent ) 

lag in seconds between signal for control and its actual motion 
time , seconds 

reel part of root of stability equation 

time for oscillation to damp to one-half its -amplitude, seconds 
control motion as a function of time 


T period of oscillation, seconds 


DTClERMKtVriON OF CONDITIONS FO? NEURAL 


The method of determining the conditio ns for neutral stability 
is Illustrated in figure 1. The calculated phase angle of the control 
motion - nd the calcula ted ratio of the amplitude of control motion to 


trplane motion are plotted against angular frequency as shown by the 
valid -lino curves. The upper lashed curve is a plot of the experimental 
ratio of the amplitude of control motion to autopilot motion against 


an yular frequency 


lower three dashed curves are three possible 


experimental phase-angle curves for the automatic pilot. The 
intersection of the experimental mid calculated control -amplitude 
curves establishes the approximate frequency of the airplane with the 
autopilot in operation. If, as in -the case of the intermediate 


experimental phase-angle curve , the intersection of the experimental 


arid calculated phase -angle curves is at the same frequency as the 
intersection of the control-emplitu&e -ratio curves, the <airplane may 
be neutrally stable and may be expected to oscillate continuously at 
this frequency. It is, however, more usual that the intersection of 
the experimental and calculated phase-angle curves will not be at 
the same frequency as the intersection of the controi-amplitude-ratio 
curves. If the phase -rngle curves, as in one case shown , intersect 
at a higher frequency than the control-amplitude* ratio curves , the 
aircraft will be stable. If, as in the remaining case, the inter- 
section of the phase -angle curves is at a lower frequency than the 
control-amplitude -ratio curve s , the aircraft will be unstable. 

Because of the nonlinear characteristics of the control! system, it 
is generally necessary to make the experiments for different ampli- 
tudes . With a dead spot (insensitivity to small deviations ) there 


will probably be some amplitude below which the system will be unstable. 
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Calculct ed freque ncy -re soonce curves f or tie a i rcraf 1m 
application of the method to the case of aileron control of 
craft or a missile Independently 5- tali 11 zed at out all three 
the equation of motion fer determining the control movement 


In the 
an air- 
axes , 




\ 

-°) 


0 



O 


The calculated steady- state sola: 
a sinusoidal forcing function of 


0 = 0 . 


max 


sin 


(art 


0. 


(St 


Ion of the p,± r 
rail t amplitude 
areuce 1.) The 


and 0 are obtained ever the desired range of 
by trie substitution of ice for D in the enuu 


■aft in response to 

o = sin art- is 

■.allies of 0. 

'naax 

gular frequencies 
on 


0) 


5 

0 


obG 


D _+ % 
Q ' l b 


This substitution 

is equivalent t 

'O specify iru' 

motion find results 

s in the express 

lion A + :LB 

obtained -i- = 

\jj\r + and 

a _ -1 B 

U — ucu: — 

ft 


vji unloinped sinusoidal 
from which can be 

The male 0 may 


denote either a phase lag or lead; depending upon its quadrant. 

If 0 is in the third or fourth quadrant; the control lags behind 
the displacement; but if 0 is in the first or second quadrant, the 
control leads the motion of the airplane. The ratio of the amplitudes 

and } is —•« and mry be termed the control- amplitude ratio X 
rmax 

of the control system, that is, the ratio of maximum control deflection 
to maximum disp lac orient in bank. A plot of K and 0 against cd 
shows the combination of control-- amplitude ratio and phase lag or lead 
necessary to maintain fixed amplitude oscillations at any given 
frequency. These results are the calculated frequency- response curves 
due to a sinusoidal motion of the aircraft. 


Determination of equi valen t, si ne wave for the a utomatic -pi lot 
response . - The experimental frequency-response curves are obtained 
by oscillating the automatic pilot sinusoidally at various amplitudes 
through the desired range of frequencies. The control i3 assumed to 
oscillate at the same frequency as the automatic pilot but because of 
the physical characteristics of the autopilot the control motion may 
differ widely fro m tx true sine wave and may chow /.arbitrary phase, 
amplitude, or wave -form relations (fig. 2). It is necessary, therefore, 
to determine an equivalent sine -vs ve response for any arbitrary control 
motion. 
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In order to determine the equivalent nine wove for an arbitrary 
control motion, the f olievir ■» relations arc acsuried: 


momentum of the airplane 
the came interval . 


a- 

control f 

a 1 

harmonic . 

3.0 

tu-al cent: 

; 7i fr' 

ivlvalont 

qu 

ui the clr 

0 

actual co: 


Trie work done “by a north nimonic force 5(t) of frequency co upon 
a harmonic motion sin cot is proportional to 

* r r . . 

B-j / c ( t ) COG evb It 


whei'o the period of the oscillation is T = — > and B n is the coef - 

ficient of the component coo crb. (f'iou reference 2?.) This con- 
p orient of the control notion that is out of phase with the air- 
craft motion is the only harmonic of the Fourier series re presenting 
the forcing function b( t ) which contributes to the work done on the 
aircraft . The angular impulse is obtained hy integrating the curve of 
control deflection against tine ever a half cycle . This component of 
t 3 control motion in phase with the sinusoidal motion ox the aircraft, 
obtained from the second relation^ is 

I T/2 

I ‘ 

At = ;•■ j S(t) dt 
1 u 0 

tho coefficient of the component sin arc 


wuere 


.to u: vL 

or zero net ! 


The con- 


dJ tion of sero net impulse over a series of cycles may be met bv 

t "> 

adjusting the reference axis for o(t) so that ■ &(t) dt = 0. 

v ' 0 

Tho control motion may then bo expressed as the sum of the in-phase 
ar.i out -of -phase components 


A, sin ct + 3, cos art 


or 


I O o . . 

JA, 1 +3,^ sin (tot + 0) 
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Tho control *nrmlltudo ratio li is oqir 
phase lag of the system d is t an* 1 ^ 

ratio and the phase in g or leaf? are cl e leru? nod freon records taken 
of the os ci Heir lone mi plotted agaiuc i, co. In contrast to the 
calculated frequency -response curves wh ich jn.oive tho aerodyimiic 
and roans characteristics of the airorafi,, those orperimen tally , 
detemined curves will he functions of tho dead spots and various 
types of lag found in tlie con orol system. In gonerel, the behavior 
of the automatic pilot will be nonlinear and hence a family of curves 
shoving different phase *rxu ariplitude /eolations fo:c different ampli- 
tudes of disturbance will be obtained. 

Comparison of tho c a lculi b ed fre oi enr^ -re on once cur ves of the 
airc raft and the oxpc.x , inen ’ 1 'a l f r e cucnc~ - r espQn n a curves of th e 
automatic pilot .- The two sets of fre ncenoy-resaoime curves show, on 
the one hand, tho values of It and *9 necessary fer hunting at a 
given frequency and, on the other hand, the actwl values of K and 0 
obtained experimentally at tills frcnvcucy. In ardor to determine from 
these curves whether the aircraft will hunt In fii-ht, the following 
conditions must bo satisfied: 

(1) At a given frequency and amplitude the experimental values 
of If and 0 must agree with the calc* luted vluos. 

(2) The 
one at which 
condition) . 

Tie first erudition indicates that the con irol- amplitude ratio 
and phase lag or lead obtained ns a result of nil types of lag In 
the control system must agree with the combine t/ion of K and 0 
necessary for him ting to exist. The second condition i3 essential 
to prevent instability if the aircraft is displaced to amplitudes 
larger than the one at which it will hunt. Tho aircraft is stable 
at these larger • mp 1 i iudes if, at the frequency for which the 
ca.lcu3.ated and experimental control -amplitude rati os are equal, the 
calculated value of the phase lag required, for hunting is greater 
than the experimental value. In other wards, the calculated value 
of 0 Is a critical value of the lag necessary to cause the -aircraft 
to hunt. If experimental values of 0 are loss than the critical 
value, the aircraft motion is damped , whereas Instability occurs if 
the experimental value of 0 exceeds the calculated critical value 
(fif!- i). 

Illustrative case . - The equation of notion in bank of a small 
experimental aircraft tested in the Langley 7* by 10-foot tunnel 


motion must be stable for amplitudes 'larger then the 
the airplane will nunc, (as dotorminod from the first 


: o o 

to dlf- -i. py and tho 

» .1. _L 

The c on 1 . col -ampli tude 
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\T5.3 estinated as 


(O.OOOSVpD 2 -f- 0.002 


= - 0.0 


Solving the equation, for 6/0 and sutstitutixig 
the expression 


I/O 


'or D give 


Tlie resultant value £ 


-7 = 0 .000726b - 0.009?ao:L 

55 

;* K and 9 , that is. 


ho calculated frequency- 
solid lines in figure 3* 


re i -pons g curves for the aircraft, are shown a 
These curves shew that for snal.1 values of K 
steady oscillation is low and the motion will not he sustained unless 
the phase lag is large, la 1C increases, the frequency of the steady 
oscillation increases hut the phase lag recurred decreases. It is 
important to note that an automatic pilot with a constant time lag l 
would he unstable at high angular 1 frequencies since the relation 
between angular frequency, phase lag, and time lag is 0 (radians) = coZ , 

Tire experimental f re quency-resv once curves were ob tained by 
oscillating the automatic pilot at amplitudes of 10° and 20° through 
the range of desired angular frequencies . Trie pha.se lag end control- 
amplitude ra.tio for* the two amplitudes were determined from records 
similar to figure 2 and are plotted as a function of co in figure 3 
for two valves of control amplitude ratio K. For this particular 
automatic pilot, the c ontro 1 - empl i 1;ude ratio was independent of 
amplitude whereas the phase lag varied with amplitude. The results 
in figure 3 indicate that, for each control -amplitude ratio, the 
experimental values of 9 pure greater than the calculated phase lag, 


and hence the aircr. 


would be unstable. Unpublished results from 


wind-tunnel tests indicated that the motion was unstable, as predicted 
from the curves of figure 3 • 


In an effort to make the aircraft stable, the parameters of the 
automatic pilot were modified and additional wind-tunnel tests were 
performed. The conditions selected for these wind-tunnel tests were, 
however, different from these conditions for which the experimental 
frequency-response curves wore obtained and hence no direct prediction 
of the aircraft stability could be made . The results of the wind- 
tunnel tests with the modified automatic pilot Indicated a steady 


oscillation, and records taken of the tests showed that 


value £ 


o** K and 0 agreed very closely with the combination of K and 9 
determined from the calculated frequency- response curves . The circled 
test points of figure h show the combination of K and 9 for the 
cocos in which steady oscillations occurred in the roll tests. 
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. Ui XU 


aircraft 


FOR BAILED C301 hi ATIONS 


The previous calculations of the f re queue;; 
were based on the assiimption that the sinusoid'': 
aircraft is neutrally damped. It is often tool 
determine the pe.rfo:ma::ce of an automatic cor it:, 
to cause the notion of the airernf 1 to d^-j.ap .r-.t 
rate. Although no satisfactory analysi 3 of thi 
given, a qualitative indication of the rate of 
in a given case nay to obtained by cojrrparing th 
amplitude of the control to the phase and ornplL 
required to enforce a given rate of dumping. 


• no spouse curve 3 
L motion of the 
vb-i-c, however, to 
device required 
sufficiently rapid 
- problem has been 
I*. ruling to be expected 
j measured phase and 
■vie calculated to be 


Strictly speaking, the assumed exponen 02 al damping of the motion 
would require an exponential decrease in the response cf tno auto- 
pilot at decreasing amplitudes. In general, such a linear response 
c anno t be e xp e c 'to d : nd 1 J 45 nc e th e me eh o& v i 1 1. re 0 '■ 1 i re c are f ul j u&gmerv 
in its application. 


The equation of damped motion in b anlr for the illustrative case 
may be written by adding a real part 

1 

; 0.000245 (-4 + in)) 2 + 0.0024: 

i 

where is given a value as detenu 


P = 0.00092a'.) 2 - 0. 00092 :5;i 2 + 0.009264 + i( - 0.009204 4 - 0.00l352cm) 


H 9 o 

the i 

maginary root 

' (-M + 

iv) 7 

= -0.26458 

ned "by 

oho d 

e sired rate o; 

.f i lK a 

ripli a 

ude in T^/g : 

>/0 feiv 

■ec th 

3 enur es si on 


The frequency -response curves shown in figure 
values cf A varying from 0.173 to 8*31* The 
?’atic curves are only plotted for M- equal to 
are the limiting curves for the values of h 1 
indicates that the control- ampli tude ratio is ?• 
whereas the phase lag decreases as ^ increase* 
with a given control- amplitude ratio, therefore 
oscillation increases ns the phase lag is redac- 
ts to damp one -half in less than 1/7*22 second, 


wore calculated for 
control- amplitude - 
zone and 5 since they 
rvesti pa ted . Figure 5 
imoet independent of M> 
s . For a control system 
, tho damping of the 
cd. If the oscillation 
the control motion 


must lead 


mo 


ion. 
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In order to predict quantitatively the stability of the motion 
of an aircraft which damps exponentially, the experimental frequency- 
response cxtrves would have to he obtained for the condition where the 
forced oscillation of the automatic pilot also damps exponentially. 

CONCLUDING REMARKS 


A method for predicting the stahility ±n roll of automatically 
controlled aircraft hy a comparison of calculated frequency -response 
curves for the aircraft and experimentally determined frequency - 
response curves for the automatic pilot is presented. The method is 
expected to he usefi.il as a means of establishing the specifications 
of the performance required of the automatic control device for 
pilotless aircraft designed as missiles. 

Langley Memorial Aeronautical Laboratory 

National Advisory Cor out tee for Aeronautics 
Langley Field, Va. 
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Fig. 1 
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C: 


Angular frequency t cj , rad /a ns J sec 
F/yu re /.- Frequency-response cun/es for stable , unstable , 
neutrally damped motion of the aircraft . 
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Fig. 3 



Angular frequency , (V , radians / sec 

F, gure 3 .— A compart son of the calculated and experimental 
frequency-response curves of a small aircraft . 
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Fig. 5 


y ' 0/4 vj dpnj.iidujD-/oj4uoj 
^ rv 



fop ' 0' 3/6' UO 39D0CJ 





